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Abstract 

We construct a gravitational dual of a Josephson junction. Calculations on the 
gravity side reproduce the standard relation between the current across the junction 
and the phase difference of the condensate. We also study the dependence of the 
maximum current on the temperature and size of the junction and reproduce familiar 
results. 



1 Introduction 



One of the most remarkable results to emerge from string theory is gauge/gravity duality 
[TJ [2J E] the equivalence between a theory of gravity and a nongravitational theory. This 
equivalence is called holographic since the nongravitational system lives in a lower dimen- 
sional space than the one with gravity. In recent years this duality has been successfully 
applied to condensed matter systems [HE]. In particular, a gravitational dual of a super- 
conductor has been found (6j [Tj. In this paper, we extend this construction to obtain a 
gravitational dual of a Josephson junction. 

A Josephson junction [8] consists of two superconductors separated by a weak link. If 
the link supports a phase difference 7 between the condensate of the two superconductors, 
there will be a current across the junction given by 

J = JmaxSin7. (1.1) 

This current exists with no applied voltage. There are several types of Josephson junctions 
depending on the nature of this link. It can be an insulator (SIS junctions), a normal 
conductor (SNS junctions), or even a very narrow superconductor (e.g. bridges and point 
contacts). 

We will construct a gravitational dual of an SNS Josephson junction and reproduce 



(1.1). We will also show that J max decreases exponentially with the size of the weak link 
as expected, and study its temperature dependence. Although our main motivation is to 
extend the range of condensed matter phenomena that can be described gravitationally, 
refinements of this approach may provide new insights into Josephson junctions made with 
high T c superconductors. 



2 The Model 



We begin with the action 

S = I d A x \f^g 



m 



(2.1) 



where F = dA and D = V — iqA. We will restrict our analysis to the probe approximation. 
That is, we rescale ip — i/)/q, A — A/q and take q — > 00, keeping ip and A fixed. In this 
limit, the backreaction of these fields on the metric can be ignored. Therefore, we can fix 
the metric background to be the AdS planar Schwarzschild black hole: 
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where ro is the black hole horizon radius and L the AdS length scale. The temperature of 
this black hole is given by 

T =T I T-2- ( 2 -3) 
4ttL 2 

This simple model has been shown to provide a dual description of a superconductor [7] . 
As a consequence of a scaling symmetry, the critical temperature T c of the superconductor 
is proportional to the chemical potential //. We will allow the chemical potential to have 
spatial dependence^ fi(x) so that a fixed temperature T is above the critical temperature 
in a narrow gap but below the critical temperature elsewhere. In this way, we model two 
superconductors separated by a normal phase. 

Without loss of generality, we will create our junction along the x direction. Our fields 
will now be functions of r and x. A Josephson junction requires a phase difference, so we 
must include a phase in our scalar field. Therefore, we consider solutions of the form 

4> = \i\)\e iLp , A = A t dt + A r dr + A x dx , (2.4) 

where tp, A t , A r , and A x are all real functions of r and x. Instead of A, we wish to work 
with the gauge- invariant fields M — A — dip. 
The equations of motion are 

d 2 .M t + 4*% M t + - d rMt ~ = , (2.5b) 

d 2 x M r - d r d x M x - 2r 2 \ip\ 2 M r = , (2.5c) 

d 2 r M x - d r d x M r + f - (d r M x - d x M r ) - 2 ^M X = , (2.5d) 

1 If M \ f f 2\ 

d r M r + ^jd x M x + — ^M r d r \^\ + -jjdM) + (j + -j M r = • (2.5e) 

Note that the phase ip no longer appears in the equations of motion. The first four equations 
(2.5a 2.5d) are second-order dynamical equations while the last one (2.5e) is a first-order 
constraint equation coming from the conservation of the source in Maxwell's equations. 

Now we discuss the boundary conditions, beginning with r = oo. At this point, we will 
work in units in which L — 1. For simplicity, we choose^Jm 2 = —2. With this choice of mass, 
the scalar has the following asymptotic form near r = oo: 



1 / f 2\ 1 / /If 2 M 2 

d 2 M + -^d'M + lL + -)d r \^\ + -l^-fM 2 -^-m 2 )\^\=0, (2.5a) 



^ (1) (x) ^ 



2), 



X 



/y>2 \ ty 3 



W = — + + • ( 2 - 6 ) 



1 For other inhomogeneous holographic superconductors, see [Hllinj: and see jTTj for a related solution. 
2 Our analysis can be extended to any mass above the Breitenlbhner-Freedman bound, m 2 > —9/4. 
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The ip^ and ip^ terms are both normalizable, so we have a choice of boundary conditions. 
Here, we only consider the case ip^ 1 ' = 0. Via gauge/gravity duality, ip^ then gives the 
expectation value of a dimension two operator in the boundary field theory 



(O) = ^ 2) . (2.7) 



We interpret this operator as the superconducting condensate. 
The Maxwell fields take the asymptotic form 



■V, -/,(,■) + o(±) . (2i 



M r = O ( \ ) , (2.9) 



M x = u(x) + - + 0[\) . (2.10) 



The quantities [i(x), p(x), u(x) and J are interpreted in the boundary field theory as the 
chemical potential, charge density, superfluid velocity, and current, respectively. By solving 
the equations of motion near r = oo, it can be shown that J must be a constant. This 



condition is imposed by the constraint (2.5e), which in turn enforces current conservation. 



Our boundary conditions are determined by choosing J and /t(x). 

At the horizon r = Tq, regularity requires M t = 0. This and the equations of motion 
place boundary conditions on the remaining functions. 

At x = ±oo, we require that the functions approach the homogeneous (x-independent) 
solution. We solve the ordinary differential equations numerically, imposing regularity at 
the horizon and choosing J and //(oo) = /i(— oo) at r = oo. With non-zero^] J, there are 
typically two solutions. We choose the solution with the lower free energy, which is also the 
solution that has a larger value of the condensate [T3] . 

The gauge invariant phase difference is 7 = Aip — f A x where the integral is across the 
gap. Since the edges of our gap will not be completely sharp (for numerical reasons), it is 



more convenient to use (2.10) and set 

7 = -/ dx [v(x) - K±oo)] ■ ( 2A1 ) 



The second term is needed as a regulator since the homogeneous superconductor at x — ±00 
has a nonzero current and hence a nonzero superfluid velocity. Note that the value of this 
term at +00 and —00 is the same. For numerical convenience, and because the phase 



5 For a discussion of holographic superconductors with nonzero current, see [T2l[T31[Tl] 
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difference 7 between the two superconductors is a non-local quantity, we will reproduce (1.1) 



by first picking J and then computing 7 using (2.11) 



The equations of motion (2.5) are invariant under the following scaling symmetry: 



r — > ar 



(t,x,y) (t,x,y)/a , M r -> M r /a , (M t , J14) ->■ a(M t , Mj) . (2.12) 



For homogeneous superconductors, the scale invariant temperature is T/fx, so changing \i 
is equivalent to changing T. In the inhomogeneous case, T must still be constant, even 
though fj, = fi(x). The parameter a in (2.12) must be constant, so there is no symmetry 



which exchanges an inhomogeneous chemical potential with an inhomogeneous temperature. 
The scale invariant temperature is T/T c where T c is the critical temperature of the junction, 
which is identical to the critical temperature of a homogenous superconductor with zero 
current and \i = //(oo) = //(— 00). T c is proportional to //(oo) and given by 



T c « 0.0588 //(oo) 



(2.13) 



When performing numerics, we will use (2.12) to set tq = 1, and adjust T/T c by varying 
/i(oo). 

We would like to choose a profile for //(#) such that our system resembles an SNS 
Josephson junction. Therefore, we choose fi(x) to be small and approximately constant 
for x G (—4, 4) and then rise rapidly to //(oo) outside this gap. Inside the gap, the effective 
critical temperature is 

T « 0.0588 //(0) , (2.14) 

For temperatures T < T < T c , the gap is a normal conductor and the region outside the 
gap is superconducting, so we indeed have an SNS Josephson junction. For temperatures 
smaller than T , our material is everywhere in the superconducting phase, and we have an 
S-S'-S junction. For T > T c , our material is entirely in the normal phase. 
A profile that fits this description is given by 



^(x) = // c 



2tanh(£ 



tanh 



x + 



a 



tanh 



x — 



a 



(2.15) 



where /too is the chemical potential at x = ±00, and i is the width of our junction)^} The 
quantities a and e control the steepness and depth of our profile, respectively. For t 3> a, 
n(x) is quite flat inside the junction (see Fig. 1). For this profile, T = eT c . 



4 We have tried a different profile which is less steep and consequently easier for numerics. This profile 



was also able to reproduce (1.1 1 
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3 Numerics and Results 



We can now attempt to solve the coupled system of nonlinear partial differential equations 



(2.5) numerically for any value of J and any parameters in (2.15). In order to proceed, it is 
convenient to impose boundary conditions at x = and x = +oo rather than at x — ±oo. 
By using the equations of motion and their perturbations, one can prove that our boundary 
conditions require that M r be odd and \if)\, M t , and M x be even. This in turn imposes, at 
x — 0, Dirichlet boundary conditions on M r and Neumann boundary conditions on the other 
functions. We can also compactify both the radial coordinate r and boundary coordinate x 
using the change of variables z — 1 — tq/t and x = tanhfx/ (4er)]. 

Our numerical method relies on a standard relaxation procedure, combined with spectral 
methods defined on a Chebyshev grid. Due to the exponential convergence characteristic of 
spectral methods on such a grid, we only require a few points in our spatial discretization 
scheme. In all the plots in this manuscript, we use 41 points along the x direction and 
25 along z. We have varied the number of points, and found good agreement with the 
aforementioned exponential convergence. 



After solving the equations, the phase difference can be obtained through (2.11). We can 
study the dependence on the size of our junction and on the temperature by varying t and 
/Zoo, respectively. As an example, we show on the left panel of Fig. [I] a typical result of our 



numerical code for M t , showing the profile (2.15) imposed on the boundary. 



A standard property of Josephson junctions is (1.1). We only consider phase differences 
that lie in the interval (— 7r/2,7r/2) since, within a cycle, the other phase differences corre- 
spond to dynamically unstable junctions [16J. By feeding into our code successive values 



of J and computing 7 using (2.11), we constructed the graph on the right panel of Fig. [T] 



The black solid line represents the best fit of our numerical data to fll.l ). The agreement is 



striking, and predicts a maximum current across the junction of J ma , x /T^ fa 1.408. We have 
used other profiles for n(x) to model the junction, and the resulting curve shares the same 
features as the right panel of Fig. [TJ showing that (1.1) is a robust signature of holographic 



SNS Josephson junctions. Two natural questions, which we address next, is how J max varies 
with the junction width i and the temperature T/T c . 

The dependence of J max on £ is shown on the left panel of Fig. [2j Once again, this curve 
is in good agreement with condensed matter physics [T3], which for SNS junctions, predicts 
an exponential decay with growing i in the maximum current: 

Jmax/T c 2 = A e~l. (3.1) 

This universal law is only expected to hold as long as £ <C I. Furthermore, a similar behavior 



5 



1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

7 

Figure 1: In the first graph, we represent Mt as a function of x and z for J/T^ « 0.0081, showing 
the boundary profile fi(x) imposed at z = 1. The second graph shows how J/T^ varies 
as a function of 7. The solid line is the best fit sine curve. In both plots we use /i^ = 6, 
e = 0.6, and a = 0.5. The first plot has £ = 6 and the second has £ = 3. 




is expected of the condensate within the barrier at zero current: 



(O) x=0 ,j=o/T* = A 1 e~5i. 



(3.2) 



Note that the value of the coherence length £ should be the same in (3.1) and (3.2). We 
have developed an independent code that studies the limiting case of zero current, in which 
and M t are the only non-zero variables. We show, as an inset plot of the left panel of 
Fig. |2| how (C) x=0 j =0 varies with I. We fit both sets of data to (3.1) and (3.2) and find 



{£,A } « {1.17,17.96} and {i,A x } « {1.26,33.52} for J max /T 6 2 and (O) x=0 ,j =0 /T?, respec- 
tively. Note that these two sets of numbers were computed using different numerical codes, 
each depending on a different number of variables, and each having their own numerical 
error. Furthermore, the normal and superconducting phases in our junction are not cleanly 
separated - our numerics prevent us from choosing a profile for fj,(x) that is too steep. We 
believe this justifies the 7% disagreement between the two estimates of £. 

We now turn our attention to the dependence of J max on the temperature. Since there 
is not a simple form for the expected behavior of J ma x(7 n ), we do not include fits to our plot 
on the right side of Fig. [2] 

There are two distinct regimes for how J max varies with T, corresponding to the two 
temperature scales (2.13) and (2.14) set by the profile (2.15). For temperatures smaller than 
T , our material is everywhere in the superconducting phase and we have an S-S'-S junction. 
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Figure 2: In the first graph, we represent J m ax and (O) x= o (inset plot) as functions of I for /ioo = 6.0. 

The second graph shows how J ma , x /T^ varies as a function of T/T c , for I = 3. In all the 
plots, we use e = 0.6 and a = 0.5. The black solid lines are the best fit exponential 
curves. 



In this regime, we find that our system deviates far from (1.1). Indeed, this behavior is 
expected of S-S'-S junctions [T7], but the deviation creates difficulties in finding J max so we 
do not include that region in our plot. In Fig. |2j we have chosen e = 0.6, so this region 
corresponds to T/T c < 0.6. 

There is a complementary temperature regime near T c , corresponding to the critical tem- 
perature at which the entire junction is in the normal phase. Although this is an interesting 
regime, it is difficult to probe numerically due to large scale separations between the con- 
densate size (O) and the charge density p. However, we observe that as T approaches T c , 
Jmax approaches zero, as is expected. 



4 Discussion 

We have seen that the simplest model of a holographic superconductor can also be used 
to construct a holographic Josephson junction by choosing /i(x) appropriately. Since the 
chemical potential is not constant, one might wonder why we find stationary solutions. 
After all, there is a nonzero electric field E(x) = fi'(x) on the boundary near the edges 
of the gap which might be expected to make p time dependent. However this effect can 
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be balanced by diffusion [18] leading to the stationary configuration we findj^] We expect 
that holographic Josephson junctions with constant /i can be constructed by modifying the 
boundary condition for the scalar field in a narrow gap along the lines discussed in [19] . 

We have restricted attention to the DC Josephson effect. The AC effect says that if a 
voltage difference V is applied across the junction, then the solution is indeed time dependent 
and the phase difference satisfies 

d 1 /dt = qV (4.1) 

so one gets an alternating current. The added time dependence makes it harder to construct 
the appropriate gravitational solutions. However, there is a good reason to expect that the 
AC Josephson effect will also be reproduced on the gravity side: it can be obtained by 
doing a bulk gauge transformation on half of the solution. Adding a voltage difference V 
corresponds to adding a constant V to A t , but this is gauge equivalent to changing the phase 



of our charged scalar ip by qVt. This is precisely (4.1). 

Another natural extension of our model is to add a magnetic field. This should cause the 
phase to increase linearly in a spatial direction along the junction. The dual gravitational 
solutions now depend on three variables. We leave a discussion of their properties for future 
work. 

Our analysis has been restricted to the probe approximation. To go beyond the probe 
limit and study the effects of backreaction, one must solve the full dynamical equations, 
including Einstein's equations. The first step in this direction was done in [14], where 
the homogeneous superconductor with nonzero current was constructed. The backreacted 
junction would require solving a set of nine coupled, nonlinear partial differential equations. 

Our model can be extended to accommodate a junction whose weak link is a narrow 
superconductor (a Dayem bridge). One would need a profile fi(x,y) that breaks supercon- 
ductivity in the regions x G (— §, §), \y\ > |. As in the case with a magnetic field, the 
gravitational solutions will now depend on three variables. We expect these junctions to 
support a larger J max - 
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